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I.  IjB  TB0D0CTI01 


There  are  aany  solutions  to  the  problem  of  predicting 
the  location  of  an  underwater  energy  source.  One  common 
solution  is  the  use  cf  a  passive  hydrophone  which  detects 
the  pressure  waves  radiating  from  the  source.  The  hydro* 
phone  sensor  is  assumed  to  be  omnidirectional  and  therefore 
incapable  of  estimating  direction.  To  provide  direction¬ 
ality  a  series  cf  sensors  are  placed  in  a  row  to  form  a 
passive  linear  array. 

■  A  familiar  methcd  of  determining  directionality  is 
time-domain  beamforming.  In  this  principle,  it  is  assumed 
that  the  source  is  far  enough  away  so  that  the  pressure  wave 
appears  to  be  a  plane  wave  when  viewed  at  the  site  of  the 
receiving  array. 

Thus  a  set  of  time  delays  are  calculated  for  any  direc¬ 
tion  cf  signal  arrival,  which,  when  applied  to  the  receiver 
outputs  causes  them  to  be  in  phase  and  to  reinforce  when 
summed.  The  resultant  angular  response  to  signals  arriving 
from  other  than  the  nominated  direction  is  then  a  function 
cf  the  array  geometry,  relative  to  the  signal  wavelength, 
and  any  weighting  factors  which  have  been  applied  to  the 
receiver  outputs.  The  effect  is  to  generate  a  main 
receiving  beam  in  the  desired  direction,  with  a  series  of 
undesired  subsidiary  side  lobes  whose  magnitude  can  be 
controlled  to  some  extent  by  the  choice  of  a  suitable  array 
geometry  and  the  use  cf  amplitude  weightings  on  the  receiver 
outputs. 

In  order  to  determine  location,  three  or  more  such 
arrays  separated  by  a  known  amount  may  be  used. 

This  study  is  concerned  with  a  single  linear  vertical 
passive  array  and  the  determination  of  the  depth  of  a  single 


under nater  source.  The  analysis  is  based  on  the  following 

assumptions: 

•  The  underwater  source  is  emitting  continuously  and  at  a 
monochrcaatic  freguency. 

•  Both  the  source  and  receiving  hydrophones  are  stationary 
in  space  causing  a  constant  range. 

•  The  range  is  sufficiently  long  so  that  the  channel  is 
filled  with  R-R  (ref  rac  ted-refracted)  rays. 

•  There  is  no  distortion  introduced  in  the  propagating 
medium  so  that  the  signals  received  at  each  sensor  are 
identical  except  for  constant  delays. 

•  The  source  signal  and  noise  are  independent  and 
stationary  gaussian  random  processes. 

•  The  speed  of  soujd  profile  is  triangular  and  symmetric 
with  the  deep  sound  channel  axis  at  1000  meters.  The 
velocity  gradient  is  -0.017  meters/meter /sec  above  1000 
meters  and  +0.017  meters/meter/sec  below  1000  meters. 
This  profile  gives  a  speed  of  sound  at  the  surface  of 
1500  meters/sec. 

•  The  speed  of  sound  profile  is  constant  in  the  horizontal 
plane. 

•  Only  R-R  rays  are  considered.  All  other  rays  have 
sufficient  loss  that  their  effect  is  negligible. 

As  opposed  to  conventional  time-domain  beamforming,  this 
study  makes  no  assumption  of  planar  wave  fronts  at  the 
receiver  site.  Therefore  the  time  delays  applied  to  each 
receiver  will  not.  in  general,  be  a  linear  function  of 
depth. 

Since  it  is  desired  to  determine  whether  or  not  there  is 
a  source  present  at  a  specific  depth  the  result  will  be  a 
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binary  decision.  A  "1"  will  indicane  signal  source 
present;  a  "0"  will  indicate  signal  source  non  present.  For 
the  constant  range  there  will  be  “N"  tesn  depths  investi¬ 
gated  for  the  signal  source.  The  number  of  hydrophones  in 
the  vertical  array  will  be  ML". 

Per  a  single  source  at  a  given  depth,  the  travel  time  is 
calculated  from  the  depth  to  each  hydrophone.  This  travel 
time  is  converted  into  a  phase  delay  for  each  hydrophone  so 
that  after  summation  from  all  hydrophones  a  maximum  output 
is  achieved.  This  output  is  then  passed  through  a  squaring 
device,  an  integrator,  and  a  threshold  and  flip  flop  device 
to  give  a  w  1"  binary  output.  If  the  signal  source  is  at  a 
different  depth  and  the  same  previous  phase  delays  are  used 
for  each  hydrophone  then  the  output  will  be  somewhat  less 
than  the  previous  maximum.  The  difference  in  depth  required 
to  achieve  a  "0"  binary  output  is  the  depth  resolution  of 
the  system. 

The  travel  times  for  each  hydrophone  are  calculated  for 
each  cf  the  "S"  source  depths  to  be  considered.  HN"  will 
ordinarily  be  much  greater  than  "L"  so  that  the  system  will 
be  overdetermined.  An  over  determined  system  is  one  in  which 
there  are  more  equations  than  unknowns.  The  objective  then 
is  to  calculate  the  phase  angle  and  amplitude  weight  for 
each  hydrophone  so  that  a  determination  can  be  made  indi¬ 
cating  the  presence  or  absence  of  a  signal  source  at  a  given 
dept  h . 

The  method  used  to  calculate  the  phase  and  amplitude 
weights  is  the  linear  minimum  variance  estimation  technique, 
linear  minimum  variance  estimators  are  optimum  when  compared 
with  all  ether  estimators  for  gaussian  problems.  The  method 
is  directly  applicable  to  overdetermined  systems. 

The  output  of  the  summer  is  calculated  using  the  linear 
minimum  variance  aiplitude  and  phase  angles  assuming  a 
source  at  one  of  the  "NM  source  depths  and  no  source  at  the 


11 


ethers.  The  calculation  is  repeared  for  each  of  the  depths. 
The  result,  when  plotted  against  source  depth,  will  be 
referred  to  as  the  ’'beam  pattern"  of  ♦he  array  in  this 
report.  (Although  si  mi  la  r ,  it  should  not  be  interpreted  as 
the  angular  response  of  an  array  as  in  the  conventional 
definition  of  a  beam  pattern.  The  conventional  definition 
loses  ouch  of  its  utility  when  the  wavefronts  are  not 
planar.)  Ideally  the  beam  pattern  will  b6  maximum  at  the 
desired  depth  and  very  small  at  ail  other  depths  so  that  the 
binary  "1"  decision  will  be  made  for  a  source  at  the  desired 
depth,  and  a  "0"  decision  for  sources  at  all  others.  This 
beam  pattern  is  compared  with  the  depth  beam  pattern  of  the 
conventional  time-domain  beamformer  mentioned  above.  The 
purpose  of  this  iVesis  is  to  determine,  as  an  initial  inves¬ 
tigation,  whether  the  linear  minimum  variance  estimation 
technique,  when  applied  to  a  linear  vertical  array,  is 
useful  in  depth  discrimination  at  long  ranges  in  a  ' SOFAR' 
type  sound  channel. 


IX.  GBBBRAL  THBOBY 


1.  BIT  ACOUSTICS 


The  propagation  c£  sound  in  an  elastic  median  can  be 
described  mathematically  by  solutions  of  the  wave  equation 
using  the  appropriate  boundary  and  medium  conditions  for  a 
particular  problem.  The  wave  equation  relating  the  acoustic 
pressure  'p*  to  the  coordinates  'x'/yS'z',  and  the  time 
*t* ,  may  he  written  as 


,2  2  .2  .2  ,2 

d  p  m  c  dp  +  dp  +  dp 

dt2  dx2  dy2  dz2 


(2.1) 


where  'c*  is  a  quantity  that  has  the  general  significance  of 
sound  velocity  and  may  vary  with  the  coordinates. 

One  may  approximate  the  solution  of  the  wave  equation 
using  ray  theory:  its  body  of  results  and  conclusions  is 
called  ray  acoustics. 

Officer  [Hef.  1]  describes  the  ray  solution  as  a 
complete  solution  tc  any  particular  propagation  problem 
within  the  validity  of  the  approximation  of  the  Eikonal 
equation  to  the  wave  equation.  For  these  approximations  to 
be  valid  neither  the  amplitude  of  the  wave  nor  the  speed  of 
sound  can  change  appreciably  in  distances  comparable  to  a 
wavelength. 

ThU3  the  path  of  a  ray  through  a  medium  in  which  the 
speed  of  sound  varies  with  depth  can  be  calculated  by  the 
application  of  Snell's  law 


cose/c  *  1/cQ  *  a  constant  for  any  one  ray  (2.2) 
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where  '€*  is  the  angle  of  depression  made  with  -hs 
horizontal  at  a  depth  where  the  spaed  of  sound  is  *  c* ,  and 
*co*  is  tha  speed  at  a  depth  (real  or  extrapolated)  where 
the  ray  would  become  horizontal. 

In  a  medium  in  which  the  velocity  of  sound  changes 
linearly  with  depth  the  sound  rays  can  be  shown  to  be  arcs 
of  circles,  that  is,  to  have  a  constant  radius  of  curvature. 
Kinsler  et  al.  [Ref.  2]  give  a  simple  and  heuristic  demons* 
tation  of  the  circularity  of  rays  in  a  medium  with  a  linear 
sound  speed  gradient  ’g*.  The  center  of  the  circle  which 
creates  the  arc  lies  at  a  depth  where  the  sound  speed  extra* 
polates  to  zero.  To  understand  this,  consider  a  portion  of 
a  ray  path  with  a  local  radius  of  curvature  *R»,  as  illus¬ 
trated  in  Figure  2.1.  Since  the  gradient  'g'  for  this  case 
is 

g*Ac/nz  *  (^-Cj^  /  (d2~d^)  *  (c2“C1)  /B  (cose1*cos©2)  (2.3) 

where  •  a  c*  is  the  change  in  sound  speed  and  z'  is  the 
change  in  dapth.  It  can  be  seen  that  the  radius  of  curva¬ 
ture  is  given  by 

R  ■  - Cq/ g  *  -c/(g  ccsfl)  (2.4) 

The  ray  path  is  therefore  a  circle  when  'g*  is  constant 
because  'B'  is  then  constant.  The  center  of  curvature  of  a 
circle  lies  at  the  depth  where  0  is  90  degrees,  which  corre¬ 
sponds  to  c*0.  For  the  situation  in  Figure  2.1  the  speed 
gradient  is  negative  so  that  fB*  is  positive.  If  the  speed 
gradient  were  positive  *R*  would  be  negative,  and  the  path 
would  curve  upward. 

Once  the  radius  cf  curvature  of  each  segment  of  a  path 
is  known  the  actual  path  can  be  traced  graphically  or 
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computed.  If  the  initial  angle  of  depression  of  a  ray  is 
01,  then  by  referring  to  the  geometry  of  Pigure  2.1  the 
changes  in  both  range  and  depth  are 

hr*  ^  (sin  (91)  -sin  (02)  )/<g  cos(01))  (2.5) 


Az*  (cos (02)  -cos  (01) )/( g  cos  (01))  (2.6) 

The  sign  convention  for  these  equations  is:  downward,  to  ~he 
right,  and  depressicn  angles  below  the  horizontal  axis 
positive. 

The  symaetric  triangular  sound  speed  profile  assumed  in 
the  introduction  is  siailar  to  speed  profiles  encountered  in 
the  deep  sound  channel,  sometimes  called  the  SOFAS  channel. 
The  velocity  minimum  which  occurs  at  the  axis  of  the  sound 
channel  causes  the  sea  to  act  like  a  kind  of  lens;  above  and 
below  the  minimum,  the  velocity  gradient  continually  bends 
the  sound  rays  toward  the  depth  of  minimum  veloci-y.  A 
portion  of  the  power  radiated  by  a  source  in  the  deep  sound 
channel  accordingly  remains  within  the  channel  and  encoun¬ 
ters  no  acoustic  losses  by  reflection  from  the  surface  and 
bottom.  These  cays  are  called  B-B  (r efracted-refr acted)  and 
since  they  have  very  low  transmission  loss,  very  long  ranges 
can  be  obtained  from  a  source  of  moderate  acoustic  power 
output.  Thus  an  energy  source  at  a  specific  depth  will 
propagate  energy  in  all  directions  but  only  the  direction 
which  is  toward  the  receiving  array  and  for  which  the  rays 
arm  B-B  is  of  interest.  To  determine  the  range  of  depres¬ 
sion  angles  which  will  yield  B-B  rays  Snell's  law  is  used  to 


9sax  *  arc  cos  (cl/ (cl- (g  d) ) ) 


(2.7) 


where  '8max'  is  in  radians  and  'cl'  is  the  speed  of  sound  at 
the  scurcs  depth  »d'. 

An  example  is  given  in  Figure  2.2  of  a  single  ray  trace 
propagating  in  the  SOFAS  channel  to  show  how  aquations  2.5 
and  2.6  are  used  in  deteraining  range  and  depth.  The  ray 
path  is  broken  up  into  arcs  of  circles  as  shown  in  the 
figure,  and  then  by  paying  close  attention  to  the  previously 
defined  sign  conventions,  the  change  in  range  and  depth  is 
found.  Being  aore  specific,  for  arc  1,  91  and  62  are  both 

positive;  fcr  arc  2,  61  is  positive  and  92  is  zero;  for  arc 
3,  01  is  zero  and  02  is  negative;  and  for  arc  4,  01  is 

negative  and  9  2  is  zero.  By  keeping  a  running  total  of  all 
the  depth  and  range  changes  it  is  possible  to  deter  nine  the 
total  horizontal  distance  travelled  and  the  depth  at  that 
distance. 

Fcr  the  spaed  of  sound  profile  assumed  in  the  introduc¬ 
tion  a  computer  generated  ray  plot  is  shown  in  Pigure  2.3 
for  a  source  depth  of  300  meters.  As  each  ray  propagates 
out  from  the  source  the  triangular  channel  becomes  filled 
with  sound.  .If  a  receiving  hydrophone  is  placed  a  great 
distance  away,  a  number  of  refracted  propagation  paths  will 
exist,  each  having  a  different  travel  time  and  crossing  the 
channel  axis  at  different  intervals.  The  path  with  the 
greatest  excursion  from  the  axis  will  have  the  shortest 
travel  time. 

Officer  fHef.  1]  shows  that  the  travel  time  't'  of  a 
ray,  which  is  an  arc  of  a  circle,  is  given  by 

t  2  1,  c2  d0 

g  ^  cosS"  (2.8) 

9 1 

and  the  travel  time  for  each  arc  is 
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(2.9) 


t 


1  loge 
g 


tan((7r/4)  +  (92/2)) 
tan((ir/4T+(S1/2)) 


Equation  2.9,  applied  using  the  sane  convention  as  equations 
2.5  and  2.6,  determines  the  total  travel  time  of  a  ray  in 
the  deep  sound  channel. 


B.  1BHAY  HODEL 

The  receiving  linear  vertical  array  is  presumed  to 
consist  of  •!'  hydrophones  as  shown  in  Figure  2.4.  It  is 
assumed  that  the  source  is  emitting  energy  at  a  constant 
frequency,  'f',  and  amplitude,  'A1,  regardless  of  the  depth. 
The  source  signal  at  the  source  is  Aexp(j27r  ft).  The 
inherent  received  signal  at  the  first  hydrophone  is 


x  1  (t)  «  A'exp  <  j2*f  (t-tl))  (2.10) 

where  'tl*  is  the  travel  time  from  the  energy  source  to  the 

first  hydrophone  and  *A* '  is  the  amplitude  of  the  signal  at 

* 

the  range  cf  the  array.  After  passage  through  the  amplitude 
weight  '  al*  and  a  phase  delay  of  »ti*,  the  signal  on  the 
first  hydrophode  at  the  input  to  the  summer  is 


yl(t)=A'alxl(t-Ti  )=Af  alexpC  j  2irf  (t-tl-Ti ) )  (2.11) 

A  time  delay,  for  monochromatic  signals  corresponds  to  a 
phase  shift 


e  *  2 ir  ft 


(2.  12) 


where  '9'  is  the  phase  shift  in  radians  and  *t*  is  the  time 
delay  in  seconds.  Thus  equation  2.11  can  be  written  as 


18 


^v:^v&C:v:v:v:vC 


yl  <t,  »A'a1exp<  j<2  ft-^1-91)) 


(2. 13) 


where  '^1*2Trft1'  is  the  phase  delay  due  to  the  travel  tine 
from  the  source  to  the  first  hydrophone  and  '91*  is  the 
phase  dalay  in  the  receiver  on  the  first  hydrophone. 

Combining  all  the  hydrophones  in  the  array  in  a  summer 
gives  as  an  expression  for  the  array  output 


U 

I(t)  »  I  i*akaxp  ( j<2*ft-$k-9kn 
1c**  1 


(2.14) 


where  *^k'  represents  the  phase  delay  due  to  the  travel  time 
from  tha  source  to  the  "It  th"  hydrophone  and  1 9^ '  is  the 
phase  delay  in  the  receiver  on  tha  "k  thn  hydrophone.  If 
the  amplitude  of  the  energy  sourca  is  normalized  by  setting 
A'*1,  and  equation  2.14  is  written  in  terms  of  real  and 
imaginary  components,  we  have 


YCt)=  k=lk  cos(-V9k)+jsin(-V9k)  (2.  15, 


When  an  enargy  source  is  at  the  "q  th"  depth,  we  wish  to 
have  each  receiving  hydrophone's  phase  delay  cancel  out  the 
effect  of  the  travel  time  from  tha  source  to  it,  such  that 
*”$k~®k'  to  a  multiple  of  '2**.  This  will  put  all 

signals  into  the  sutmer  in  phase  and  thus  maximize  the 
signal  gain  for  a  source  at  tha  "q  th"  depth.  From  equation 
2.15, 


I^Ca^  cos^”^kq”®kq^  =  L  (source  present  at  q)  (2.16) 


u 

1^(3^  Sin^“^]cq“9kq^  5  9  (source  present  at  q) 


(2.  17) 


i^wiNVyv-v.v:-.:  .>>>::  -• »» 


Note  that  the  first  subscript  on  the  phase  angle  indicates 
the  receiving  hydrophene  and  the  second  subscript  indicates 
the  depth  of  the  source.  Thus  *^q  '  would  indicate  the 
phase  shift  relating  to  the  travel  time  free  the  "g  th"  test 
depth  to  the  Mk  th"  hydrophone  in  the  receiving  array. 

It  is  desirable  for  *  Y(t)'  to  be  a  minimum  value  for 
sources  at  other  than  the  depth  being  investigated.  Thus 
for  each  of  the  other  '  H-1*  depths  »Y  (t)  1  is  set  to  zero. 
This  gives  '8-1*  equations  for  the  real  terms  of  'Y(t)'  set 
to  zero 


£  a^cos  a0  (source  absent;  nt/q)  (2.18) 

k=l 

and  •  8-1  •  equations  for  the  imaginary  terms  of  'Y(t)'  set  to 
zero 


iA<is 


in  (*tknf  9kq’ 


(source  absent;  m/q) 


(2.  19) 


By  using  elementary  trigoncmetric  identities,  equation 
2.16  (real  terns  with  source  prasent  at  "q  thM  depth) 
becomes 


k-1  *kq  [C°S  ^k<l,COS  (flkq)-sin  (fkq)  sin  (©kq>]  *1-  (2.  20) 

Equation  2. 18  (real  terms  with  source  absent  for  each  of  the 
ether  *8-1  *  depths)  becomes 

Ji*1Cq[co»^jan)co»(9kq)",in(l>taa,8in<9kq)]  s0  ®sX’2’ *  *  (2.21) 

Equation  2.17  (imaginary  terms  with  source  present  at  the  "q 
th"  depth)  becomes 


20 


y  -a,  (Sin  (4  )cos(6  )+ccs(<k  )sin(9  ))=0  (2.22) 

k“x  kq  Tkq  kq  Tkq  kq 

and  equation  2.19  (imaginary  terms  with  source  absent  for 
each  of  the  other  *  H-1*  depths)  becomes 


u  p 

J1-akq[8in^tan)cos(8kq)+G°s^)an>»in(9k<l)]  =0  *=1,2, ...N;  m*q  (2.23) 

Thus,  there  are  a  total  of  '2N*  equations  with  '  2L' 
unknowns. 

In  order  to  simplify,  we  put  these  * 2N •  equations  into 
matrix  fcrm.  Arbitrarily  the  real  terms  are  made  the  first 
•H *  equations  and  the  imaginary  terms  the  second  ' S'  equa¬ 
tions.  The  first  real  and  first  imaginary  equation  is  at 
the  lowest  (shallowest)  source  depth  and  equations  increase 
in  order  after  that  until  the  last  real  and  last  imaginary 
equation  correspond  tc  a  source  at  the  deepest  depth.  The 
resultant  matrix  equation  becomes: 


0 

COS  4 1  COS$71  ...  COS*L1 

>-sin*11  -sin*21  ...  -sin*L1 

alqC°89lq 

0 

• 

cos$12  cos$22  cos^L2 

■•sin$22  •  •  • 

a2qCOs92q 

• 

0 

L 

0 

co»*lq  cos*2q  * ‘ ’  C°8*Lq 

-sin^lq  -sin*2q  ...  -sin*Lq 

0 

costlf,  cos*2N  ...  cos*LN 

-sin^1N  -sin^2jj  ...  -sin^^ 

aLqC089Lq 

0 

-sin*^  -sin*21  .-sin*^ 

-cos*^  -cosi>21  ...  -cos*L1 

alq8in9lq 

0 

-sin412  -sin$22  .-sin^j 

-cos$12  -cos^22  ...  -cos^L2 

a2q8in92q 

0 

-sin^m  -sin^2M  . -sin$LM 

-cos*1N  -cos*2fJ  ...  -cos«»LM 

aLq8in9Lq 

Simplifying  further,  equation  2.24  becomes: 
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2.25 


where  *  *jg_*  ^  *£.' #  and  '.sj  represent  ths  appropriate 

sub a a trices. 

Then,  fcy  noting  that  the  multiplication  of  each  element 
of  a  column  by  the  same  nonzero  constant  doesn't  affect  the 
solution,  equation  2.25  becomes: 


0  -JrrtTl  H 

finally  letting  Z,  and  0  represent  the  matrices  in  equa¬ 
tion  2.26,  we  obtain 

Z.  *  10  (2.  27) 


where  a  matrix  is  denoted  by  a  capitalized  underlined 
letter. 

In  summary,  £  is  the  desired  response  of  the  vertical 
array  to  the  'S'  source  test  depths.  A.  represents  known 
travel  times  from  each  of  the  ' N •  source  depths  to  each  of 
the  'L'  receiving  hydrophones.  is  unknown.  It  is  the 
phase  and  amplitude  weighting  which  must  be  applied  tc  the 
vertical  array  in  crder  to  realize  £,.  ^.contains  *2L' 
unknowns. 


Equation  2.27  represents  *2 N '  equations.  Since  this 
system  of  equations  is  over  determined  (H>L)  an  exact  solu¬ 
tion  does  not  exist.  In  order  to  make  the  best  estimate  of 
g  for  the. desired  response,  the  linear  minimum  variance 
estiaaticn  technique  is  used. 

C.  LIIEAB  BUI  HUH  Y1BXAHCE  HETBOD 

Equation  2.27  represents  a  noise  free  environaent.  If 
noise  were  present  it  would  become 

Z  *  A0+n_  (2.  28) 

with  'n.'  a  "2H"  eleaent  coluan  vector.  This  represents  the 
noise  at  each  source  depth.  is  a  linear  function  of  *6 1 

and  is  called  the  oberserva ticn.  4'  is  a  w2Nx2L"  modula¬ 

tion  or  observation  matrix  which  is  known,  and  *£•  is  the 
"2L"  element  random  parameter  vector  which  is  to  be  esti¬ 
mated.  Assume  the  first  and  second  moments  of  •©•  and  *ja* 
are  given  by 

E(£)  58  iie  Y ar (9)  *  Y0  (2.29) 

and 

E(n)  *  0  Var(n)  *  Vn  (2.30) 

where  *E'  represents  expectation  or  first  moment  and  ’Var' 
represents  the  covariance  aatrix.  It  is  assumed  that  the 
paraaeter  '6*  and  the  noise  '.a.'  are  uncorrelated. 

A  restriction  iaposed  is  that  the  estimate  must  be  a 
weighted  linear  combination  of  the  observations: 


(2.31) 


iL  - 

where  ,/v'  indicates  estimate.  The  objective  is  tc  select 
'b/  and  *B_*  in  order  tc  minimize  the  error  variance.  Such 
an  estimator  is  the  linear  minimum  variance  estimator;  it  is 
the  test,  in  the  sense  cf  minimum-error-variance  linear 
estimators. 

Another  restriction  is  that  the  estimator  be  unbiased: 
in  otherwcrds  it  is  required  that  the  expected  value  of  the 

A 

estimator  •©  *  is  equal  to  the  expected  value  of  the  param- 
L 

eter  Thus, 


E(Sl)  =t>BE(z)  =  E<0)  =  Hq  (2.32) 

yielding 


b=ii9-EAUQ  (2.33) 

Substituting  this  result  in  equation  2.31  gives  for  the 
unbiased  linear  estimator 


®L3^0+£<-~-H.e)  (2.34) 

Hote  that  since  the  estimator  is  unbiased,  the  estimation 
error  '0  =6-6  *  is  zero  mean.  The  next  step  is  to  select 

C«  Li 

•J.'  in  order  to  minimize  the  error  variance.  However,  this 
optimization  problem  is  ill-defined  because  the  error  vari¬ 
ance  is  a  matrix.  Therefore  in  order  to  introduce  a  scalar 
goodness  measure  the  sum  of  the  variances  of  each  component 
of  is  minimized.  This  is  the  sum  of  the  main  diagonal 
terms  of  the  covariance  matrix  and  is  defined  as  the  trace 
cf  the  matrix. 


tr  (Var(e  ) )  =  I  Var((9r)  ) 
~L  n=l  *  n 


(2.35) 


«hare  'tr'  indicates  trace.  '£'  is  then  selected  to 
minimize  the  trace  of  the  error  variance,  cr 


ain  tr  (7ar(0_))  =  ain  tr  (E  (0  ei,)  ) 
B  “E  B  L  h 


(2.36) 


where  ,T '  indicates  the  transpose.  The  following  problem  is 
then  obtained  by  substituting  equation  2.34  into  equation 
2.36: 

tain  tr  ( Var(a„) )  =min  tr(E(  ^9-0^ -3(2-^ )  ( I-UgBCZ-AUg  )T) )  (2.37) 

3  3 

It  is  well  known  (  Hef,  3]  that  equation  2.37  is  minimized 
when 


cov  (9 , Z)-BVar  (Z)  =0 


(2.  38) 


where  'Ccv(0,Z)'  is  the  covariance  matrix  of  the  unknown 
parameters  and  the  observations.  Denoting  the  optimum 
filter  by  'J,*' ,  then  if 


B*=Cov  (0,Z)  (Var  (Z)  )  -1 


(2.39) 


a  minimum  is  achieved  for  the  sum  of  the  squares  of  the 


errors. 


Using  equation  2.28  for  '£' ,  the  covariance  of  and 
becomes 


Co  v  (0 ,  Z)  =  Cov  (0 , A0+n)  =Vg  &  1 


(2.40) 


since  and  are  uncorrelated.  The  variance  of  'Z_ '  is 
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Substituting  these  into  equation  2.39  gives 


* 

B 


*  VeAT(M6AT+Vn)-1 


and  the  linear  minimum  variance  estimator  is 


(2.  42) 


— LMV 


Ee AT  ( AVg  AT + )  “ 1  (  2  -  AUq  ) 


(2.43) 


By  utilizing  a  matrix  inversion  lemma  [Bef.  3]  equation 
2.43  becomes 


iLMV  *  <AT^1A+^1>'1<ATv;1Z*v;1ue)  (2.44) 

Th9  advantage  of  this  equation  over  aquation  2.43  is  the 
size  of  the  matrix  to  be  inverted.  In  equation  2.4  3  the 
matrix  has  dimensionality  *  2N'  while  in  equation  2.44  its 
dimenionality  is  only  '  2L*  .  Thus  the  advantages  of  the 
linear  variance  estimator  are  the  ease  with  which  they  are 
derived,  the  mathematical  tractability  of  the  linear  form, 
and  the  minimum  amount  of  stochastic  information  required 
for  development.  &n  interesting  characteristic  is  that  the 
linear  minimum  variance  estimate  is  the  orthogonal  projec¬ 
tion  of  '&*  onto  the  space  spanned  by  the  observation  1  z*  . 
Eecause  of  these  factors  this  estimator  is  a  popular  form 
for  estimating  unknowns  in  cverdetermined  equations. 

For  this  thesis  it  is  assumed  that  the  noise  samples  are 
uncorrelated  and  identically  distributed  so  that: 


2 

a  I 


(2.45) 
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v-v:v>^v_ 


.  •  -  e  ,  -  ,  •  .  •  .  "  *  »  »  •  ,  •  v*  .  * 


1 


No  previous  knowledge  is  assumed  about  '8'  .  This  implies  an 
infinite  variance  matrix  which  is  represented  as: 


and  Ug 


(2.46) 


The  linear  minimum  variance  estimate  given  by  equation  2.44 
is  then 


®T, mv  (-l]  _1iTi  (2  •  47) 


By  determining  '  6  •  the 

"”LMV 

found  for  a  signal  source  on 


!i 

V 

alqC0S9lq 

42qCOa92q 

jjt 

3 

Lq  UJ 

9L*1 

4lq3in9lq 

>- 

aLq9in9M 

Upon  solving,  this  equation 


phase  and  amplitude  weights  are 
the  '  q  th'  depth.  Recall  that 

(2.  48) 

ives  for  the  phase  delay 


A  A  A 

9  aarctan  (9.  V9  )  where  m»i  to  L 
mq  L+m  m 


and  amplitude  weight 


(2.49) 


a  *e  /cos  (0  )  where  m=1  to  L  (2.50) 
mq  m  mq 

Mhen  these  amplitude  weights  and  phase  delays  are  applied  to 
the  vertical  linear  array  a  resulting  beam  pattern  is  formed 
which  in  the  absence  of  noise  is: 
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The  resulting  beaa  pattern  'Z  '  can  then  b®  compared  with 
the  desired  beam  pattern  ' ZJ ,  as  wall  as  with  a  conventional 
beaa  pattern  1 Z 1 1  obtained  using  linear  phase  shifts  across 
the  array  aperture  selected  to  "steer”  the  array  to  the 
dominant  arrival  angle  for  the  selected  source  depth. 
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figure  2.4  lrray  Uodel 


HI.  1IPBBI HBHTAL  PBOCEDOBB 


A.  BASIC  ASSOHPTIOBS 

The  speed  of  sound  profile  given  in  the  Introduction  was 
used  to  test  the  techinigue,  with  'L*,  the  number  of  source 
depths  chosen  as  20.  The  source  depths  are  at  220  meters 
and  every  20  meters  thereafter  to  600  meters  inclusive. 
'N ' ,  the  number  of  hydrophones  used  in  the  vertical  array 
was  chosen  as  5.  The  hydrophones  are  placed  at  depths  of 
100  ,  200  ,  300  ,  4  0  0,  and  500  meters.  Thus,  there  are  40 

equations,  two  for  each  source  depth,  and  10  unknowns,  twc 
for  each  hydrophene.  A  range  of  200  kilometers  (km)  assures 
that  the  deep  sound  channel  is  filled  with  sound  over  the 
aperture  of  the  vertical  linear  array.  A  frequency  of  100 
Bz  provides  good  resolution  in  the  beam  pattern  without 
introducing  alias  mainlobes  at  the  selected  range  across  the 
depths  of  invest igaticn. 

B.  •A*  IATBH  CALC  OI  ATI  OH 

Since  the  gradient  of  the  sound  velocity,  *g*  is 
constant  in  the  area  of  the  source  depths,  and  the  speed  of 
sound  at  the  surface  is  known,  equation  2.3  is  used  to  solve 
for  'c2',  the  speed  of  sound  at  the  source  depth.  Equation 
2.7  is  used  for  each  source  depth  to  calculate  the  maximum 
initial  depression  argle  which  yields  R-R  rays.  In  deter¬ 
mining  travel  time,  cnly  depression  angles  from  each  source 
which  are  positive  (downward)  and  yield  R-R  rays  are  used. 

Beginning  with  the  first  source  depth  of  220  meters,  an 
initial  depression  angle  of  0.0  degrees  is  selected.  The 
ray  path  is  then  calculated  using  aquations  2.5  and  2.6  and 
broken  into  a  series  of  arcs  as  in  figure  2.2.  Equation  2.9 
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is  used  to  determine  the  travel  tiae  for  each  arc  ir.  the 
saae  aancer.  Each  arc's  horizontal  range ,  travel  time,  and 
depth  are  saaaed.  Bhen  the  suamed  horizontal  range  reaches 
200  ka  the  saaaation  process  ends.  The  travel  tiae  and 
depth  of  the  ray  at  this  horizontal  range  is  then  known. 
The  saae  procedure  is  repeated  for  an  initial  depression 
angle  of  0.1  degrees  and  every  increment  of  0.1  degrees 
thereafter  until  the  zaxiau  ■  depression  from  equation  2.7  is 
reached.  The  final  ray  path  is  at  this  maximum  depression 
angle. 

The  saae  procedure  is  repeated  for  each  of  the  other  19 
source  depths. 

Thus,  for  each  initial  angle  froa  each  source  depth 
there  is  a  ray  which  has  a  travel  time  and  a  depth  when  it 
reaches  the  horizontal  range  of  200  ka.  Since  it  is  the 
profile  of  the  sound  pressure  wave  which  impinges  on  the 
vertical  array  which  is  of  importance,  a  constant  can  be 
subtracted  froa  these  calculated  travel  times.  This 
constant  is  selected  to  be  the  travel  tiae  for  the  source 
depth  of  220  meters  which  has  an  initial  depression  angle  of 
0  degrees.  It  is  subtracted  froa  each  of  the  travel  times 
making  the  resultant  travel  times  relative  with  respect  to 
the  ray  which  has  a  0  degree  depression  angle  froa  the  220 
meter  depth.  The  program  and  its  listing  which  calculates 
the  relative  travel  times  and  the  depths  of  these  rays  at 
the  horizontal  range  of  200  ka  is  given  in  Appendix  A. 

A  plot  of  the  relative  travel  times  versus  depth  for  the 
220  meter  source  depth  is  shown  in  Figure  3.1.  Figure  3.2 
displays  the  plot  for  the  3  80  meter  source  depth.  Negative 
relative  travel  tiaes  in  the  plots  indicate  that  the  overall 
travel  tiae  is  less  than  the  reference.  These  rays  arrive 
at  the  200  ka  horizontal  distance  before  the  reference  ray. 

Since  the  receiving  hydrophones  are  at  set  vertical 
positions  (100,  200,  300,  400,  and  500  aeters)  ,  an 
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interpolation  is  done  to  determine  relative  travel  times  to 
then  from  each  scarce  depth.  The  interpolation  program  and 
its  listing  is  in  Appendix  B.  sometimes  more  than  one  R-R 
ray  travels  from  the  source  depth  to  a  hydrophone.  when 
this  occurs,  the  ray  which  arrives  first  is  used  in  the 
calculation  of  relative  travel  time  to  that  hydrophone. 

Bguation  2.12  determines  the  phase  shift  relating  to  the 
relative  travel  times.  The  'JJ  matrix  is  formed  by  taking 
the  appropriate  sine  and  cosine  values  as  in  equation  2.24. 
The  ’.A.*  matrix  is  *40  by  10'. 

C.  BATBIX 

Referring  to  equation  2.26,  the  matrix  is  a  *40  by 

1'  column  vector.  It  is  the  desired  beam  pattern.  The 
bottom  20  rows  give  the  imaginary  terms  and  are  set  to  zero. 
The  top  20  rows  represent  the  value  of  the  real  terms  at 
each  source  depth.  Therefore  each  of  the  top  20  rows  is  set 
to  zero  except  for  the  row  containing  the  source.  It  is  set 
to  1 .  Fcr  example,  if  the  source  is  at  220  meters  then  only 
the  tcp  rcw  is  set  tc  1 .  If  the  source  is  at  380  meters 
then  only  the  ninth  row  is  set  to  1. 

0.  BES0LTI1G  BE  AH  P1TEEBI 

1  •  Osinq  the  Linear  Minimum  Vari ance  Het hod 

A 

•eTU„  '  is  calculated  using  equation  2.47.  The 

—LMV  a 

resulting  beam  pattern  *jZ'  is  calculated  using  equation 
2.51.  The  program  which  calculates  the  'A'  matrix,  uses  it 

A  ""  A 

in  determining  'Slmv*'  and  then  calculates  is  given  in 
Appendix  C.  The  program  listing  is  also  included. 


2 .  Osjjig  Linear  Ehase  Shifts 


The  conventional  beaa  pattern  is  determined  by  using 

A 

equation  2.51  where  '5.'  is  calculated  by  approximating  the 

plot  of  relative  travel  tine  vs.  depth  by  a  straight  line  at 

the  receiving  hydrophone  depths.  For  example  Figure  3.3 

represents  this  plot  for  the  380  meter  source  depth.  The 

straight  line  is  determined  by  a  least  squares  linear 

regression  which  minimizes  the  sum  of  the  squares  of  the 

deviations  of  the  actual  data  points  from  the  straight  line 

of  best  fit.  Note  that  only  data  points  which  are  on  the 

dominant  curve  are  used  in  calculating  the  straight  lire. 

From  the  straight  line,  relative  travel  times  tc  the 

receiving  hydrophones  are  calculated.  The  relative  travel 

tines  for  the  3  80  meter  source  depth  are  given  in  Table  I. 

They  correspond  to  a  plane  wave  arrival  angle  of  3.73 

degrees.  is  determined  by  converting  these  relative 

travel  times  to  phase  delays  using  equation  2.12  and  then 

taking  the  appropriate  sine  and  cosine  values  of  these  phase 

delays  as  in  equation  2.24.  The  amplitude  weights  are 

initially  assumed  to  be  unity. 

&  second  method  for  obtaining  the  conventional  beam 

pattern  is  calculated  by  the  same  procedure  except  the 

amplitude  weights  which  are  determined  by  equation  2.50,  the 

•0  '  amplitude  amplitude  weights,  are  applied  to  each 

LMV 

hydrophone. 


T1BL1  I 

Relative  Travel  Tiees  For  380  Beter  source  Depth 


Hydrophone  Depth 
100  outers 
200  outers 
300  aeters 
400  aeters 
500  aeters 


Relative  Travel  Time 
-0.07565509  sec. 
-0.07131599  sec. 
-0.06697690  sec. 
—0.06263780  sec. 
-0.05829871  sec. 


Pig  ace  3.1  Relative  Travel  Tiae  vs.  Depth  (220  aster  source) 


Belative  Travel  Tiae  vs-  Depth  (300  aster  source) 


Figure  3.3  Straight  Line  Approx,  of  Rel.  Trav.  Ti«e  vs.  Depth  (source  at  380  ■) 


IT.  HBSOLTS 


A.  EXACT  SOLOTIOH 

An  exact  solution  is  derived  for  the  bean  pattern  if  the 
number  of  receiving  hydrophones  equal  the  number  of  source 
depths.  For  example,  when  the  5  receivers  are  used  to 
discriminate  between  5  source  depths  (220,  240,  260,  280, 

and  300  meters)  there  are  10  equations  with  10  unknowns. 
Figure  4.1  is  a  plot  of  the  resulting  beam  pattern  with  the 
source  at  the  shallowest  depth.  Note  that  because  of  roun¬ 
doff  errors  in  the  'IHSL*  subroutines  there  is  a  small  value 
for  the  resulting  beam  pattern  at  the  non -energy  source 
depths  under  investigation. 

B.  POOR  DEPTHS  HITH  THO  RECEIVERS 

Fcr  source  depths  at  220,  280,  300,  and  320  meters  with 

receiving  hydrophones  at  100  and  200  meters,  there  are  8 

equations  with  4  unknowns.  Pigure  4.2  is  a  plot  of  the 
resulting  team  pattern  with  a  source  at  the  220  meter  depth. 

Pigure  4.3  is  the  plot  for  the  source  at  the  280  meter 

depth . 

C.  THEHTI  DEPTHS  HITH  THO  BECEITERS 

For  all  20  source  depths  with  receiving  hydrophones  at 
100  and  200  meters,  there  are  40  equations  with  4  unknowns. 
Figure  4.4  is  a  plot  of  the  resulting  beam  pattern  for  a 
source  at  the  220  meter  depth.  Figure  4.5  is  a  plot  for  the 
source  at  the  380  meter  depth. 

In  order  to  determine  if  the  'S'  calculated  in  this  case 
is  the  best  an  alternate  method  is  devised.  Four  source 


40 


depths  (220,  240,  260  meters,  and  another  source  test  depth) 
with  the  original  source  at  the  220  meter  depth  and 
receivers  at  100  and  200  aaters  are  used.  The  beam  pattern 
is  calculated  each  tiae  with  a  different  source  test  depth 
substituted  for  the  fourth  source  depth.  The  5  best 
resulting  beaa  patterns  are  selected  along  with  the  8  source 
depths  (220,  240,  260  meters,  and  the  5  test  depths  which 
created  the  5  best  keaa  patterns)  .  Then,  using  these  3 

a* 

source  depths,  is  determined  for  the  two  receivers. 

A 

This  *QJ  is  applied  to  the  two  receivers  and  the  beam 
pattern  obtained  for  all  20  source  depths. 

The  resulting  beaa  pattern  obtained  by  this  alternate 
method  isn't  as  good  as  the  beam  pattern  obtained  by  using 

A 

all  20  source  depths  in  the  determination  of  '0'. 

D.  T8E1T1  DEPTHS  HITH  PITH  RECEIVERS 

Per  all  20  source  depths  with  all  5  receiving  hydro¬ 
phones  there  are  40  equations  with  10  unknowns.  Figure  4.6 
is  a  plot  cf  the  resulting  beaa  pattern  with  the  source  at 
the  220  meter  depth.  Figures  4.7,  4.8,  and  4.9  are  the 

plots  for  the  source  at  the  360,  380,  and  400  meter  depths 

respectively. 

E.  CGVVEITIOIAL  BEAHFORHER 

Figure  4.10  is  a  plot  of  the  beam  pattern  for  a  conven¬ 
tional  beamformer  using  linear  phase  shifts  across  the  array 
with  the  source  at  the  38  0  meter  depth  and  the  amplitude 
weights  set  to  unity.  fill  20  source  depths  and  5  receiving 
hydrophones  are  used.  Note  that  in  figure  4.10  that  ‘here 
is  less  than  1  db  discriaation  between  each  of  the  source 
depths.  Figure  4.11  is  the  plot  obtained  for  the  amplitude 
weights  set  to  values  determined  by  equation  2.50. 


V.  BilGZ  OF  250  KILCBETBRS 


The  calculations  were  repeated  for  a  range  of  2  50  ha 
using  the  sane  20  source  depths  and  5  receiving  hydrophones. 
Figures  4.12  and  4.13  represent  the  plots  of  relative  travel 
tiaes  versus  depth  fcr  the  220  and  380  aetsr  source  depths 
respectively.  Figure  4.14  represents  the  straight  line 
approximation  of  the  relative  travel  tiaes  for  the  380  aeter 
depth.  Note  that  in  this  figure  the  relative  travel  tiaes 
are  represented  by  twc  straight  lines;  the  upper  line  repre¬ 
sents  linear  phase  shifts  of  the  slower  travel  tiaes  for  the 
ccnventicnal  beaaforaer  while  the  lower  line  represents 
linear  phase  shifts  of  the  faster  travel  times.  Tables  II 
and  III  are  the  straight  line  interpolations  of  these  slower 
and  faster  travel  tiaes  which  correspond  to  arrival  angles 
of  4.04  and  -3.16  degrees  respectively.  Figures  4.15  and 
4.16  are  the  resulting  beaa  patterns  for  the  conventional 
beaafcraer  for  the  slower  travel  tiaes  using  unity  amplitude 
weights  and  linear  oiniaua  variance  amplitude  weights 
respectively.  Figures  4.17  and  4.18  are  the  bean  patterns 
for  the  faster  travel  tiaes. 

Figures  4.  19,  4.20,  4.21,  and  4.22  represent  plots  of 

the  beaa  pattern  for  the  source  at  the  220,  340,  360,  and 

380  aeter  depths  respectively. 
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TIBLI  II 

Slouet  Bay  Travel  Tines  For  250  ka  Range  and  380  a 

Source 


Hydrophone  Eepth 

Relative  Tr avel 

100  aeters 

-0.09949109  sec 

200  aeters 

-0.09478615  sec 

200  aeters 

-0.09008121  sec 

400  aeters 

-0.08537627  sec 

500  aeters 

-0.08067133  sec 

i _ i 


T1BLB 

III 

Faster  Travel  Times  For  2 50 

ka  Bange  and  380  a  Source 

Hyd^O£icfle  J3egth 

Relative  T^r avel  Tine 

100  aeters 

-0. 10 100745  sec. 

200  aeters 

-0.10468763  sec. 

300  aeters 

-0. 10836781  sec. 

400  aeters 

-0.  11204799  sec. 

500  aeters 

-0.  11572817  sec. 

Figure  4.2  Beae  Pattern  (4  depths  2  receivers,  source  at  220  ■) 


ete 


n  *20  depths 


Conventional  Bean  Pattern  (source  at  380  meters) 


-unitY  a»P 


Figure  4.12  Relative  Travel  Tiee  vs.  Depth  (range=250  k»,  source  at  220  ■) 
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Figure  4.13  Relative  Travel  Tiee  vs.  Depth  (range=250  ke,  source  at  380  ■) 


SLOWER  RAYS 


a=380  m\ 


Figure  4.15  Conventional  Bean  Pattern  (H=250  ka,  d=380  n,  slower  tines 


Figure  4.17  Conventional  Beae  Pattern  (8=250  lea,  d=380  n,  faster  tines,  a=1) 


cer  ti.6s 


Figure  4.20  Beaa  Pattern  (R-250  lea,  20  depths,  5  receivers,  source  at  340  a) 
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V.  DISCUSSION  OF  BB SUITS  AHD  BECOMHEN DATIOBS 


For  nary  of  the  cases  examined,  the  linear  minimum  vari¬ 
ance  estimation  technique  gives  high  resolution  in  the 
"depth  beam  pattern"  for  sources  at  long  range.  Figures  4.2 
and  4.3  represent  beam  patterns  for  four  depths  with  two 
receivers  where  the  main  lobe  is  at  the  desired  source  depth 
with  a  beaawidth  of  30  meters.  For  20  depths  with  2 
receivers  the  system  is  highly  overiatermined  and  the  beam 
pattern  portrayed  in  figure  4.4  has  its  main  lobe  80  meters 
from  the  source  depth.  At  the  source  depth  the  strength  of 
the  beam  pattern  is  3.6  db  down.  For  figure  4.5  the  main 
lobe  is  on  the  source  depth  with  a  beaawidth  of  2  50  meters. 

For  all  20  source  depths  with  5  hydrophone  receivers  the 
results  range  from  a  beam  pattern  having  its  main  lobe  on 
the  source  depth  (figure  4.7)  with  a  beaawidth  of  20  meters 
and  a  secondary  lobe  12.1  db  down  to  a  beam  pattern  having 
its  main  lobe  40  meters  from  the  source  depth  (figure  4.22). 
For  the  team  pattern  in  figure  4.22  the  width  of  the  main 
lobe  is  80  maters  and  the  strength  of  the  beam  pattern  at 
the  scurce  depth  is  0.8  db  down. 

When  using  the  conventional  beamformer,  the  beam  pattern 
results  with  amplitude  weights  determined  by  the  linear 
minimum  variance  estimation  method  were  superior  to  the  beam 
pattern  determined  by  using  unity  amplitude  weights. 
However  in  all  cases  the  conventional  beamfcrmer  was  signif¬ 
icantly  inferior  to  the  beam  pattern  determined  by  the 
linear  minimum  variance  estimation  technique  in  Terms  of 
depth  resolution. 

The  linear  minimum  variance  estimation  technique  could 
not  be  used  for  all  ranges  with  the  chosen  sound  speed 
profile.  For  ranges  of  187  km  and  235  km  the  relative 


travel  times  produce  an  *A_'  matrix  which  is  so  ill- 

conditioned  that  the  'IMSL1  subroutine  '  LIN V2F *  cannot 

.  T 

determine  an  accurate  inverse  of  'A, 

There  is  sufficient  evidence  from  this  initial  investi¬ 
gation  that  the  linear  minimum  variance  estimation  technique 
applied  to  a  long  linear  vertical  array  can  yield  high  reso¬ 
lution  depth  information  about  passive  sources  at  very  lcng 
ranges.  However,  further  investigation  is  needed  before  any 
field  tests  are  in  order.  Recommendations  for  further 
research  are: 

1.  The  use  of  a  more  realistic  speed  of  scund  profile, 
preferably  one  actually  characteristic  of  the  ' S  CFAR* 
channel. 

2.  The  use  of  mere  hydrophones  in  the  vertical  array. 
As  more  receivers  are  used  the  beam  pattern  should 
approach  the  desired  beam  pattern  more  closely  (a 
less  over determined  system  should  yield  smaller  total 
minimum  mean  square  error)  . 

3.  Investigation  into  causes  for  the  peak  response 
falling  at  other  than  the  desired  depth  and  techni¬ 
ques  for  correcting  this  problem  when  using  the 
linear  minimum  variance  estimation  techinique. 

4.  Alternate  assumption  for  choosing  the  ray  paths  to 
include  in  the  matrix.  For  example  cne  might 
choose  the  rays  with  the  greatest  intensity  instead 
of  the  shortest  travel  time. 

5.  The  possibility  of  estimating  range  as  well  as  the 
depth  of  a  passive  source  with  linear  minimum  vari¬ 
ance  estimation  techniques. 

In  conclusion,  the  linear  minimum  variance  estimation 
technique  of  beamforming  was  significantly  superior  to  the 
conventional  beamforner.  High  resolution  depth  information 
about  passive  sources  at  long  ranges  is  provided. 
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APPBHDII  A 

RELATIVE  TR&VBL  TIHE  CALCULATIOH 


This  program  calculates  the  relative  travel  times  from 
each  source  depth  tc  the  horizontal  range  and  the  depth  of 
the  ray  at  this  range.  To  save  paper,  only  the  two  shal¬ 
lowest  source  depths  are  listed. 
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i£PMDIX  B 

IHTEBPOLATI OH  OF  RELATIVE  TRAVEL  TIHE  CALCULATION 


This  program  does  an  interpolation  of  the  output  data 
generated  in  Appendix  A.  The  relative  travel  times  are 
interpolated  for  the  receiving  hydrophone  depths. 
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RECEIVING  HYDROPHONE  DEPTHS  ARE  SELECTED.  IMSL  KOUTINc  'ICSEVU* 
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APPENDIX  C 

RESULTING  BEAN  PATTERN  POR  CALCULATED  WEIGHTS 

This  program  uses  the  output  of  Appendix  B  to  calculate 
the  'A'  matrix.  From  this,  the  amplitude  and  phase  weights 
are  determined  by  the  linear  minumum  variance  estimation 
technique.  These  weights  are  applied  to  the  array  and  the 
beam  pattern  is  obtained. 
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